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ABSTRACT 

Dual epoch spectroscopy of the lenticular galaxy, NGC 4203, obtained with 
the Hubble Space Telescope (HST) has revealed that the double-peaked compo- 
nent of the broad Ha emission line is time variable, increasing by a factor of 2.2 in 
brightness between 1999 and 2010. Modeling the gas distribution responsible for 
the double-peaked profiles indicates that a ring is a more appropriate description 
than a disk and most likely represents the contrail of a red supergiant star that 
is being tidally disrupted at a distance of ~ 1500 AU from the central black hole. 
There is also a bright core of broad Ha line emission that is not time variable and 
identified with a large scale inflow from an outer radius ~ 1 pc. If the gas number 
density is > 10^ cm~'^, as suggested by the absence of similarly broad [O I] and 
[O III] emission lines, then the steady state inflow rate is ~ 2 x 10^^ Mq /jt 
which exceeds the inflow requirement to explain the X-ray luminosity in terms of 
radiatively inefficient accretion by a factor of ~ 6. The central AGN is unable to 
sustain ionization of the broad line region, the discrepancy is particularly acute 
in 2010 when the broad Ha emission line is dominated by the contrail of the in- 
falling supergiant star. However, ram pressure shock ionization produced by the 
interaction of the in-falling supergiant with the ambient interstellar medium may 
help alleviate the ionizing deficit by generating a mechanical source of ionization 
supplementing the photoionization provided by the AGN. 

Subject headings: galaxies: Seyfert, galaxies: individual(NGC 4203), quasars: 
emission lines 



1. Introduction 



NGC 4203 is a lenticular galaxy located at a distance of 15.14 Mpc ( Tonry et al. ||200T ). 
The galaxy is of interest because it harbors a low ionization nuclear emission-line region 
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(LINER; Heckman 1980) of type 1.9 indicating that it exhibits a broad Ha emission hne 
(Ho, Fihppenko, & Sargent 1997 Ho, Fihppenko, Sargent &; Peng 1997). High angular 



resolution observations with the Space Telescope Imaging Spectrograph (STIS) aboard the 
Hubble Space Telescope (HST) revealed that the broad Ha emission line is actually double- 
peaked ( [Shields et a'I~|[2000| ) with a FWZI of ~ 12,000 km/s. NGC 4203 was observed 
with STIS again 11 years later revealing a significantly brighter double-peaked broad Ha 
emission line. Thus, NGC 4203 joins a short but growing list of LINERs with time variable 



broad Balmer emission lines of which the best studied nearby example is NGC 1097 (Storchi- 



Bergmann, Baldwin & Wilson ||1993; Storchi-Bergmann et al., ||1995, 


1997 


but also includes 
) and NGC 3065 


M81 (Bower et al. 1996 ] 


Devereux, Ford, Tsvetanov, &; Jacoby 


2003 


(Eracleous & Halpern 


200] 


). NGC 4203 also hosts a time variable, compact core, flat- 


spectrum radio source 


Nagar, Falcke, Wilson & Ulvestad, 2002 


Ulvestad & Ho 


2001). 



Maoz et al. (2005) showed that the nucleus of NGC 4203 is also variable in the ultra-violet 



(UV). However, variability has yet to be established in X-rays (Plan et al. 2010 Younes, 
Porquet, Sabra, fc Reeves ]|2011 ). Maoz (2007) compiled the radio-to- X-ray spectral energy 
distribution of the active galactic nucleus (AGN) showing that it resembles more luminous 
Seyferts. Several attempts have been made to estimate the black hole (BH) mass using 



spatially resolved gas kinematics (Shields et al. 2000 Sarzi et al. 20011 2002 Beifiori et 



al. 2009) with mixed results. Subsequently, Lewis & Eracleous (2006) adopted the stellar 



velocity dispersion measurement of Barth, Ho, & Sargent (2002) to estimate a BH mass of 
(6 ± l)x 10"^ Mq, which, when compared to the bolometric luminosity, indicates that NGC 
4203 radiates well below 



40" 



of the Eddington luminosity limit (Lewis & Eracleous 



2006 Maoz 2007). 



Accretion disks are expected to produce broad emission lines with characteristics that. 



in principle, betray their geometries (Lewis, Eracleous, & Storchi-Bergmann 2010 Gezari, 



Halpern & Eracleous 2007). For example, disks of low inchnation are predicted to produce 



an asymmetric profile with the blue peak brighter than the red (Chen & Halpern 1989). 
As the disk becomes more edge-on, the blue and red peaks achieve a similar brightness and 
may even merge into a single peak if the disk has a large ratio of outer radius to inner 
radius as may be the case in NGC 3065 ( Eracleous fc Halpern |[2001 ). However, disk models 
invariably are unable to simultaneously explain both the broad shoulders and the bright core 



of broad line emission often seen centered at the systemic velocity of the galaxy (Down et al. 



2010; Gezari, Halpern & Eracleous 2007 Eracleous & Halpern 1994), especially in nearby 



galaxies observed with high angular resolution. For example, in NGC 4579, NGC 4450 and 



NGC 4203, Barth et al. (2001), Ho et al. (2000) and Shields et al (2000), respectively, 



attribute the central core emission to a "normal" broad component. Similarly, [Down et al. 



(2010) and Gezari, Halpern & Eracleous (2007) ascribe the broad Ha emission lines seen in 
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distant quasars and broad line radio galaxies in terms of an accretion disk and a "standard" 
or "classical" broad line region (BLR). However, labeling the single peak of core emission as 
the "normal", "standard" or "classical" BLR does not necessarily entail an explanation of 
what it is exactly. There are two possibilities; it is either an inflow or an outflow. Radiatively 
driven outflows are not possible, however, in LINERs, because they radiate well below the 
Eddington luminosity limit. Thus, inflow is the only remaining explanation for the central 
core emission seen in LINERs with double-peaked profiles. 

This paper is the third in a series that attempts to model the broad emission line 
profiles for low luminosity AGNs (LLAGNs) with known central masses and derive a size 
for the BLR. Inflows can produce the single peaked broad emission line profiles seen in M81 



and NGC 3998 (Devereux & Shearer 2007 Devereux 2011) and inflows are expected to 



be associated with accretion disks, because the former, presumably, fuels that latter. An 
objective of this paper, therefore, is to demonstrate the utility of an inflow plus disk model to 
explain double peaked broad emission line profiles with central core emission and, by way of 
illustration, apply this combined model to interpret the time- variable broad Balmer emission 
lines seen in NGC 4203. The wider context for this investigation is to better understand the 
broad line phenomenon in LLAGNs by characterizing their physical properties using archival 
HST-STIS spectroscopy. NGC 4203 is an excellent candidate to study as it is nearby, bright 
and radiates well below the Eddington luminosity limit and is therefore unable to sustain a 
radiatively driven outflow thereby simplifying the interpretation of the broad emission lines 
by eliminating one of the possible explanations for their origin. 

Since the BLR is unresolved, the spatial distribution of emission cannot be directly 
measured. However, the three dimensional gravitational field strength is known, by virtue 
of the BH mass, thus the relationship between velocity and radius may be established, given 
a kinematic model for the broad emission line gas. In this way, one can, in principle, ex- 
ploit the exquisite velocity resolution of STIS to model the broad emission line profiles 
and determine a size for the BLR that heretofore has previously only been possible using 



reverberation mapping techniques (Peterson 1993, 2001). Profile fitting complements re- 



verberation mapping as the latter yields BLR sizes for a different sample of more luminous 
AGNs. 

The layout of the paper is as follows. In Section 3, the BELs seen in NGC 4203 are 
evaluated in the context of an inflow plus disk model. Some physical properties of the BLR 
are presented in Section 4. Conclusions follow in Section 5. We begin, however, with Section 
2 and a review of the emission lines observed in the nucleus of NGC 4203. 
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NGC 4203 Emission Lines 



NGC 4203 has been observed with STIS twice using the G750M grating and once 
using the G430L grating. Table 1 summarizes the observations and the archival spectra 
are presented in Fig 1. A STIS spectrum showing the broad Ha emission hne has been 



presented previously by Shields et al. (2000, 2007) based on observations obtained under 



PID 7361. Shields et ah (2000), Bassin & Bonatto (2003), Gonzales Delgado et al. (2004) 



and Sarzi et al. (2005) have all presented the G430L spectrum obtained under PID 7361. 



However, the G750M spectrum obtained under FID 11571 is shown for the first time in Fig. 
1. 

Multiple exposures obtained using the same grating and at the same position have been 
combined using the STSDAS task ocrreject. Dithered exposures were shifted using the 
STSDAS task sshift prior to combining. Subsequently, the STSDAS task xld was used to 
perform a 7 pixel wide extraction along the slit direction and centered on the nucleus. Each 



extraction samples > 80% of the encircled energy for an unresolved point source (Froffitt et 



ah 2010). 



The G750M spectra reveal a broad Ha line with a FWHM ~ 3000 km/s that brightened 
by more than a factor of two in the almost 11 year interval since the first spectrum was 
obtained. Two narrow but resolved [S II] lines appear on the red side of the broad Ha line at 
~ 6750A and two narrow but resolved [O I] lines appear on the blue side at ~ 6350A. These 
narrow lines are virtually identical in the two datasets. Thus, the brightening appears to 
be limited to just the broad Ha + [N II] emission line complex which is remarkable because 
the more recent observations were obtained with a smaller 0.1" slitQ The G430L spectrum 
reveals a wide swath of emission lines. Collectively, the G750M and G430L spectra resolve 
the Ha, H/3 and H7 lines. A more detailed description of the spectra follows, beginning with 
the Hydrogen (H) Balmer lines. 



2.0.1. Broad Ha Emission 

The Ha emission line profile exhibits broad shoulders that are accentuated in the 2010 
spectrum compared to the one obtained in 1999. Both of the [N II] vacuum wavelength 
6549.85 and 6585.28 A emission lines and the narrow component of the vacuum wavelength 



^The continuum did change slightly, down by 15% compared to 1999 in the sense that the 2010 continuum 
is slightly redder but it is not clear whether this difference is due to the different slit sizes used for the 
observations or whether it reflects an intrinsic change in the AGN. 
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6564.61 Ha emission line can be seen in a 7 pixel extraction centered on the nucleus in both 
datasets which greatly facilitates their modeling. The STSDAS contributed task specfit 
was used to model and subtract the forbidden lines and the results are reported in Table 
2 along with the fluxes for the two [S II] vacuum wavelength A6718.29 and A6732.67 A 
emission lines and the two [O I] vacuum wavelength A6302.04 and A6365.53 A emission lines. 
The flux for the model [N II] lines and the model narrow Ha line has been chosen so that 
the difference spectrum does not show an 'over-subtraction' of the broad Ha profile. The 
procedure is valid because the observed broad Ha emission line profile is otherwise smoothly 
varying. In support of the procedure, the model parameters deduced for the [N II] lines 
are consistent with the expectation based on atomic physics that sets the wavelength of the 
fainter [N II] line relative to the brighter [N II] line, constrains the flux of the fainter [N II] 
line to be 1/3 that of the brighter [N II] line, and requires that the [N II] lines share the same 
width. Two components were required to model the brightest [N II] line in order to obtain a 
satisfactory subtraction as noted in Table 2. A single gaussian was used to model the other 
lines. The brightest [O I] line is wider than the brightest [N II] line at the ~ 90% confidence 
level. That these forbidden lines have similar widths and further that the narrow Ha line is 
fainter than the brightest [N II] line are characteristics shared by other AGNs observed with 
STIS for which these lines can be clearly seen dNoel-Storr ^tal~||2003| . A model narrow line 
spectrum is illustrated in Fig. 2 along with the difference spectrum which reveals a broad 
Ha line profile that is very similar to the H/3 line described in the next section. Figure 2 
illustrates that the consequence of subtracting the forbidden emission line model is to reveal 
a broad Ha emission line symmetric about the A6588A wavelength expected for a systemic 
velocity of 1077 km/s deduced from the peak of the brightest [N II] emission line. Thus, 
there is no apparent redshift between the broad Ha emission line and the systemic redshift 
of the host galaxy. The broad Ha emission line fiux deduced from the 1999 G750M spectra 
is ~ 40% larger than estimated previously by Shields et al. (2000). The difference can be 
attributed to the fact that Shields et al. (2000) subtracted not only the forbidden emission 
lines described above, but also another component that they call the "normal" broad line 
emission centered at the systemic velocity of the galaxy. Comparing the entire Ha emission 
line flux measured in the 2010 spectrum with that measured for the 1999 spectrum indicates 
that the broad line increased in brightness by a factor of 2.2. 



2.0.2. Broad Hf3 and H^y Emission 

The H/3 and H7 emission lines are compromised only slightly by [O HI] emission lines 
which were modeled and subtracted using the STSDAS contributed task specfit as illus- 
trated in Fig. 3 and Fig. 4. The [O HI] vacuum wavelength A4364.44 A deserved special 
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attention as it is very faint and marginally resolved. A model was employed for the [O III] 
A4364.44 line in which the wavelength and line width was fixed and the flux adjusted so 
as to not over-subtract the H7 line in the difference spectrum (Fig. 4). This yielded the 
upper limit for the [O III] A 4364.44 line reported in Table 3. Subtracting the [O III] lines 
reveals broad H/3 and H7 emission lines. Unfortunately, the analysis can not reliably be 
extended to include other Balmer lines as they are even fainter and even more confused with 
emission from other ions. Consequently, emission line fluxes are reported in Table 3 for all 
the lines that can be reliably resolved and measured in the G430L spectrum including H/3, 
H7, the blend of the vacuum wavelength AA3727.09, 3729.88 A [O II] lines, the two [O III] 
vacuum wavelength A4960.30 A and A5008.24 A emission lines, and the upper limit for the 
vacuum wavelength A4364.44 A [O III] emission line. According to the diagnostic diagram 



of [Kewley et aL| ( [20061 their Fig. 5) the [O III] A5008.24 /[O II] AA3727.09, 3729.88 and 
[O I] A6324.99/Ha emission line ratios measured with STIS in 1999 qualify NGC 4203 as a 
Seyfert and the difference between the permitted and forbidden line widths further qualifies 
NGC 4203 as a Seyfert 1. 



2.0.3. Similar Balmer Line Profiles 



Figure 5 illustrates the striking similarity between the Ha and H/3 lines in NGC 4203 
when they are normalized to their respective peak intensities and the wavelength scales 
converted to velocity using the non-relativistic Doppler equation. The similarity between the 
emission line profiles is particularly impressive considering the variety of models employed 
to subtract the superimposed forbidden lines. On the other hand, the H7 line profile is 
quite different in the sense that the line core appears to be broader and the broad shoulders 
are absent, particularly noticeable on the blue side. It is not clear whether the difference 
between the H7 profile and the other Balmer lines is a consequence of lower signal-to-noise 
data for this weak line or whether it represents a real radiative transfer effect. Nevertheless, 
the broad shoulders on the Ha and H/3 lines suggest an accretion disk, as noted previously 



by Shields et al. (2000), whereas the line core is similar to what would be expected for 



a spherically symmetric shell of gas in radial motion as noted previously by Devereux fc 



Shearer 


( 


2007 


) and 


Devereux 


( 


2011 



Section 3. 
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2.0.4- B aimer Decrements 

Despite the similarity between the broad Ha and H/3 hne profiles, the Balmer decrement, 
Ha;/H/3 = 5.0 ± 0.1 is significantly different from the Case B value, 2.75, in the sense 
that the observed value is systematically 182% higher. Interpreting this ratio in terms of 
dust extinction leads to a color excess E{B — V) ^ 0.6 mag which is inconsistent with the 



reddening estimated to the AGN X-ray emission by Plan et al. 


(2010). Such deviations 


from recombination theory have been noted for other LINERs (e.g. 


Devereux ||2011 Storchi- 


Bergmann et al., 1997; Bower et al. 1996 Filippenko & Hapern 1984) and have been 



attributed to coUisional excitation in gas of high density or time variable dust extinction. 



3. The Inflow Plus Disk Model 

As alluded to in the Introduction and in Section 2.0.3, the broad Ha (and H/3) line 
profiles can be modeled using the combination of an accretion disk and an inffovij^ since 
neither model alone is able to explain the entire broad emission line. The disk model, de- 



scribed previously by Chen & Halpern (1989), describes the emission expected from a thin 



axisymmetric, relativistic accretion disk. The model includes self-consistently all relevant 
relativistic effects, such as Doppler boosting and transverse and gravitational redshifts and 
involves five free parameters, a dimensionless inner radius, ^j, and outer radius, ^o, an in- 
clination angle measured from the disk normal to the line of sight, i, an emissivity law of 
the form e oc r"'', and a velocity dispersion for the gas, a in km/s. The inflow model has 



been described most recently by Devereux (2011). It describes a steady state spherically 



symmetric infall and invokes just two free parameters, an inner and outer radius; Tinner and 
f outer 1 respectively. Collectively, the disk plus inflow model encompass 8 free parameters if 
one includes the scaling that determines the relative contribution of the disk and inflow to 
the total line emission (see Tables 4 & 5). 

Illustrative models presented in Fig. 6 show that the broad Ha line profile observed 
in 1999 is dominated by the inflow with only a small contribution from the disk, except at 
the extremities of the line. However, the situation is quite different for the broad Ha line 
profile observed in 2010 which, according to the model, appears to be dominated by the 
accretion disk, except in the line core. The inflow model parameters are chosen to be the 
same for the 1999 and 2010 data, motivated by the fact that the narrow lines, including the 



radiatively driven outflow of gas can be ruled out for NGC 4203 because the diminutive luminosity 
generated by the AGN is simply unable to provide sufficient radiation pressure to overcome the gravitational 
force of the BH. 
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narrow component of the Ha line, did not change between 1999 and 2010. Thus, only the 
disk changed. Experimentation revealed that the 2010 Ha line profile is best modeled with a 
ratio of outer to inner disk radius, ^o/6 = 3.75, a velocity dispersion, a = 350 km/s and an 
emissivity index, q=3. A more thorough exploration of the parameter space was performed 
by computing the value of xied ^ §"^1 of 156 models per epoch spanning 26 inclination^ 
between 40° < i < 90°, and 6 values for the inner radius spanning 500 < < 3000, expressed 
in units of the gravitational radius, Vg = GM,/c^. The remaining parameters are fixed as 
listed in Table 4. 



where Oj represents the observed normalized line profile intensities, Mj, the model 
normalized line profile intensities, u is the number of degrees of freedom and S the uncertainty 
in the observed normalized line profile intensities, taken to be 4% for both epochs. The 
summation was performed over the velocity span of the broad Ha line as specified in Table 5. 
The minimum Xred computed for the 2010 data was unity indicating that the axis-symmetric 
disk plus inflow model provided an excellent representation of those data. A contour plot 
of the Xred surface, presented in Figure 7, shows that there are multiple combinations of 
disk inclination and inner radius that can minimize Xred each epoch. Thus, the model 
representations of the data presented in Fig. 6 are by no means unique. However, the fact 
that the red and blue Xmin contours do not overlap suggests that there is no single disk 
model that can simultaneously explain both the 1999 and 2010 observations. There could 
even be two disks. But, the most conservative interpretation is that there is just one disk 
for which the inclination did not change between 1999 and 2010. In this case the size of the 
disk must have changed. The arrow in Fig. 7 represents the least size evolution solution 
which implies that the disk is highly inclined, i = 82° ± 1°, which presumably aided in its 
detection, and the inner radius of the disk decreased from 2770rg ± 182rg in 1999 to 2200rg 
± 216rg in 2010, the implications of which are discussed further in Section 4.7. 



■^Inclinations smaller than 40° yielded no Xmin — 1 solutions for the 2010 observations. 




(1) 
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Discussion 



4.1. BLR Size 



That the inflow plus disk model is able to reproduce the rather complex shapes observed 
for the broad Ha emission line profiles in NGC 4203 allows the sizes of the inflow and disk 
components to be estimated. The mass distribution determines that the size of the inflow is 
large; ~ 2 pc ( 27 milli-arcsec) in diameter]^ and it terminates at the much smaller ~ 14 x 
10~^ pc (200 /iarcsec) diameter disk. 

The physical size inferred for the inflow and the disk components causes NGC 4203 to 
not conform to the correlation between BLR size and UV luminosity established for quasars 



and high luminosity AGNs using reverberation mapping (Peterson 



2001 



1993 



Kaspi et al. 



2005). An extrapolation of the BLR size - luminosity relationshippj of [Kaspi et al. (2005) 



down to the low UV luminosity estimated for the AGN in NGC 4203, predicts a radius for 
the BLR that is ~ 0.5 l.d, which is ~ 20 times smaller than inferred for the inner radius of 
the disk using line profile fitting (Section 3). Conversely, the large outer radius determined 
for the inflow in NGC 4203 using profile fitting makes it physically larger than any BLR 
measured using reverberation mapping, 4 times larger than the previous record holder; the 
quasar 3C 273 (Kaspi et al. 2005). However, as noted by Kaspi et al. (2005), the BLR 
size - luminosity correlation appears to break down for low luminosity AGNs, which, with 
L(2500 A) = 4.4 X 10^° erg/s ( |Maoz ||2007[ ), would include NGC 4203. Of course, the BLR 
size - luminosity relationship of Kaspi et al. (2005) is defined by quasars and AGNs that 



are orders of magnitude more luminous than NGC 4203. Thus, the very large discrepancy 
arising from the comparison strongly suggests that the BLR in NGC 4203 is not simply that 
of a scaled down quasar. 



4.2. Virial Black Hole Masses 

Estimating the mass of the BH in NGC 4203 using the so called 'virial method' leads to a 
value that is substantially lower than the mass estimated from the stellar velocity dispersion 



by Lewis & Eracleous (2006). For example, the formalism of Greene & Ho (2005), which 



uses the FWHM and luminosity of the broad Ha emission line, underestimates the mass of 
the BH in NGC 4203 by a factor of 348 using the broad Ha emission line observed in 1999 



^With an angular diameter of 0.02" the BLR of NGC 4203 is spatially unresolved with STIS. 
^FITEXY at 1450A. 
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and a factor of 51 using the one observed in 2010. This dichotomy is regarded as further 
evidence that the BLR in NGC 4203 is very different from the BLR in more luminous AGNs. 



4.3. Broad Line Region Ionization 



Following Maoz (2007), the ionizing continuum generated by the AGN in NGC 4203 
may be represented as a power law with an optical to X-ray spectral index a = 1.11 and 
a normalization provided by the UV continuum described in Section 4.1. Integrating the 
spectral energy distribution from 13.6 eV to 100 keV using the method described previously 
in Devereux (2011 ) predicts 1.7 x 10^^ ionizing ph/s after correcting for 0.05 mag of Galactic 
extinction only. For comparison, the broad Ha emission line flux measured in 1999 (Table 
2) which is dominated by the inflow (Fig. 6) corresponds to an Ha luminosity, L{Ha) = 
1.14 xlO^ L0, which, assuming 45% of the ionizing photons are converted into Ha photons 
(Case B recombination at a temperature of 10^ K) requires (3.2 ± 0.04) x 10^^ ionizing ph/s 
using, 



Nion = L{Ha)aB/al[hvHa 



(2) 



where a^/ = 1.16 x 10 cm'^ s ^ is the effective recombination coefficient and as = 
2.59 X 10~^^ cm^ s^^ is the total Case B recombination coefficient. 

Thus, the number of ionizing photons required to excite the broad Ha emission line 
measured in 1999 exceeds the number of ionizing photons available from the AGN by a 
factor ~ 2, independent of the gas density and the filling factor. The discrepancy becomes 
larger for the brighter Ha emission line measured in 2010. Then the number of ionizing 
photons required to excite the Ha emission line exceeds that available from the AGN by 
a factor of ~ 4. Technically, these are upper limits because the AGN is variable in the 



UV and could have brightened prior to and since 2003 when it was last measured (Maoz 



et al. 2005 ^ Nevertheless, tentatively, NGC 4203 joins M81 and NGC 4579 as examples 



of LINERs with a BLR ionizing deficit ( 


Devereux & Shearer 


2007; 


Maoz et al. 


1998 


Ho, 



Filippenko & Sargent 1996). 



^Thc ionizing deficit increases by another factor of 5 if the extinction to the Ha emission line is estimated 
from the Balmer decrement mentioned previously in Section 2.0.4. On the other hand, the anomalous Balmer 
decrement may be the consequence of high gas density as noted in Section 4.7. 
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4.4. Constraints on the Inflow Gas Density 

The size of the aperture employed to extract the spectra iUustrated in Fig. 1 encom- 
passes emission from gas up to a radial distance of 13 pc from the central AGN where the 
velocity dispersion is expected to be ~ 200 km/s and comparable to that measured for the 
[S II] lines. Thus, the two [S II] lines provide an opportunity to measure the gas density 
well beyond the outer boundary of the inflow. The nominal value for the observed [S II] 
A6742/A6757 intensity ratio = 0.68 ± 0.1, corresponding to n = 2 x 10^ cm^^, with uncer- 
tainties permitting densities in the range 1 < n (10^ cm~^) < 6. The fact that the FWHM 
of the broad Ha and H/3 emission line^is larger than that of the [O I] and [O III] emission 
lines (Table 2 & 3) means that the gas density inside the inflow must be greater than the 
critical density of the levels from which the [O I] and [O III] lines originate which is ~10^ 
cm~^. This limit is higher than the one obtained from the [S II] line ratio implying that the 
gas density increases towards the central BH. 

4.5. The Mass of Ionized Gas Required to Produce the Broad Ha Emission 

Line 

The mass of emitting gas may be deduced from the broad Ha emission line luminosity 
assuming standard (Case B) recombination theory; 

Memitting = L{Ha)mH /na'^JiUfio^ (3) 

Using an effective recombination coefficient a^^ = 1.16 x 10"^'^ cm'^ s~^, assuming a 
constant average density n > 10® cm~'^, and a luminosity L{Ha) = 1.14 x 10® Lq based 
on the broad line flux measured in 1999 (Table 2), leads to an upper limit on the mass of 
inflowing ionized ga^ Memitting < 10 M0. If only a fraction of the gas is ionized then the 
upper limit on the total (ionized + neutral) gas mass could, of course, be higher. 



^Observed in 1999 and hence dominated by the inflow. 

^The brighter hne in the 2010 spectrum approximately doubled the mass of ionized gas compared to that 
seen in 1999 if the number density is the same. Although, evidence presented in Section 4.7 suggests that 
the number density for the gas producing the double-peaked broad emission line in 2010 is much higher. 
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4.6. An Assessment of the Inflow Scenario for NGC 4203 

It is straiglit forward to calculate the volume filling factor, e, for the ionized gas pro- 
ducing the Ha emission once the dimensions of the emitting region have been established. 
For a uniform density medium occupying a spherical volume of radius r, one finds 



3L(/7„)/47rn^a^^„hz/Har' 



(4) 



Again, using an effective recombination coefficient a'^^ 



Ha 



1.16x 10 



-1 s - 
cm s 



assuming 

a constant average number density n > 10^ cm~^, and a luminosity L{Ha) = 1.14 x 10^ L0 
based on the 1999 measurement of the broad line ffux reported in Table 2, leads to an upper 
hmit on the fiUing factor, e < 1 x 10-^ for NGC 4203 if the size of the BELR, r = 1 
pc. The very low filling factor suggests that the infiow is not continuous but composed of 
many ionized, density bounded, gas filaments. Such filaments would have approximately 
the same gas density and hence the same emissivity regardless of their location with respect 
to the central AGN, they would be optically thin and hence emit isotropically. Ionized gas 
filaments embrace all the essential elements of the infiow model and are commonly seen in 



the nuclei of active galaxies ( Storchi-Bergmann 2009) including the central parsec of the 



Galactic Center (Lacy, Achtermann & Serabyn, 1991). Consequently, the infiow modeled 



here is not a Bondi fiow (Bondi 1952) because the infiowing 'particles' are discrete and do 



not constitute a continuous fiuid. 



Having established the dimensions of the emitting region and a lower limit on the gas 
density one can estimate the mass infiow rate, m, for the ionized gas observed in 1999, using 
the equation of continuity; 



m 



eAnr vnmn 



(5) 



The free-fall velocity, v, at a radius of 1 pc is determined by the mass distribution to be 
730 km/s. Setting e < 1 x 10~^ for a self-consistent gas density in the fiow of n >10^ cm~^ 
(see section 4.4), one obtains a mass infiow rate, rh ~ 2.4 x 10~^ Mq/jt. However, if only 
a fraction of the infiowing gas is ionized, then the total mass infiow rate could, of course, be 
higher. 

The 1-10 keV X-ray luminosity adopted for the AGN in NGC 4203, Li_io kev, is 7.3 x 
IQ^o gj,g g-1 ( Pian et al. [2010 ). Assuming this is powered by radiatively inefficient accretion 



leads to the following formula (Merloni, Heinz, & Di Matteo 2003) 
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10 keV 



7 X 1038M°-9^m2-3 



(6) 



where Li_io keV is in erg and M, is in solar masses. Under these circumstances 
the accretion rate required to power the observed X-ray emission, fn ~ 3.9 x 10~^ M© yr~^. 
Thus, the steady state spherically symmetric inflow inferred from the broad Ha line emission 
seen in 1999 exceeds the requirement to explain the X-ray luminosity in terms of radiatively 
inefficient accretion by a factor of ~ 6. 



The most obvious source for the inflowing material is stellar- mass loss (e.g. Ho 2009) 
with the caveat that it is only the gas for which the vector cross product 



rxp = 



(7) 



where r is the radius vector and p is the linear momentum of the gas, will be accreted by 
the BH. This is expected to limit the accretion rate to the values < 0.1 yr~^ measured 



for the inflows in M81 ( [Devereux k Shearer ||2007| ), NGC 3998 ( [Devereux |[20Tl| ) and NGC 
4203 studied here. 



4.7. An Assessment of the Accretion Disk Scenario for NGC 4203 



There are four indications that the accretion disk modeled in NGC 4203 is not, in fact, 
an accretion disk. First, the FWHM and the total broad line flux both increased between 
1999 and 2010 (Table 2) which is contrary to the expected behavior for a disk subjected to a 



brightening of the central AGN (see for example, the discussion of 3C 332 in Gezari, Halpern 



& Eracleous 2007). Second, there was likely no associated brightening of the central AGN 



as discussed further in the latter part of this section. Third, the disk is more appropriately 
described as a ring given the small ratio of outer to inner radius deduced from the model 
(Table 4). Fourth, the inner radius of the ring apparently decreased in size by ~ 20% 
between 1999 and 2010. Collectively, these observations suggest that we are witnessing a 
tidal disruption event unfolding in NGC 4203. Tidally disrupted solar-type stars are unlikely 
candidates as the inner radius of the ring (1687 AU in 1999 and 1340 AU in 2010)[^is about 
a factor of 10'^ larger than the tidal radius expected for a 1 Mq main sequence star in the 
vicinity of a 6 x 10^ M0 BH. However, the tidal radius increases for objects of lower density. 



^There are 0.609 AU/r^ in NGC 4203. 
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For example, the tidal radius, Rj at which co-rotating objects of density, p are disrupted by 
a BH of mass, M, is given by 



i?i = O.78([p0/p][M./lO^M0])^/=^ AU (8) 

Thus, the average density of an object that will be tidally disrupted at a distance similar 
to the inner radius deduced for the ring (Section 3), corresponds to ~ 10~^ pQ, comparable 
to the average density expected for a late-type supergiant star, similar to Betelgeuse for 
exampl^^ Consequently, it is quite plausible that the brightening of the double-peaked 
broad lines in NGC 4203 was caused by the ongoing tidal disruption of a red supergiant star 
which would have made ~ 1.5 orbits in the time span of the ST/S" observations (see Davidge 



2010, for evidence that red supergiant stars exist in SO galaxies). 



In the tidal disruption scenario, the average density of ionized gas responsible for the 
double-peaked broad Ha emission line is likely to be comparable to the outer atmosphere 
of Betelgeuse for which n ~ 5 x 10^^ cm~^ (Ravi et al. 2011), high enough, perhaps, to 



explain the anomalous Balmer decrement noted in Section 2.0.4. The mass of gas can be 
estimated using equation 2, an effective recombination coefficient = 1.16 x 10^^'^ cm^ 
s~^ and a luminosity L{Ha) = 1.38 x 10^ Lq, based on the difference between the broad line 
flux measured in 2010 and 1999 which should more closely represent that produced by just 
the ring in 2010. The resulting mass^^ of ionized gas corresponds to ~ 2 x 10~^ Mq and is 



presumably strewn in a ring shaped contrail marking the trajectory of the supergiant much 



as envisaged by Scoville & Norman (1995, see their Fig. 1) and Bogdanovic et al. (2004, and 



references therein). More recently Maiolino et al. (2010) describe objects shaped like comets 
orbiting and eclipsing the central X-ray source in NGC 1365. The defining characteristics 
of their "comets" - a dense head and low mass tail - are intriguingly similar to the those 
proposed here for the supergiant producing the stellar-contrail in NGC 4203. 

It was noted in section 4.3 that NGC 4203 suffers from an ionizing deficit that was 
particularly acute in 2010. The deficit could be alleviated if the AGN had brightened by a 
factor of 4 since 2003, but this is unlikely as the amplitude of the UV variability is ~ 50% at 
the shortest wavelength when it was measured by Maoz et al. (2005). Furthermore, one can 



not appeal to a boost of X-rays as no X-ray variability has been detected (Plan et al. 2010 



^"The tidal radius increases for objects of lower density and the average density of Betelgeuse could range 
from 10^^ pq to 10"^'^ pq depending on its radius which could be as large as 134 AU 

^^It is quite remarkable that the capabilities of HST are such that it can detect radiation from an ionized 
stellar contrail with a mass equivalent to 686 Moon masses in a galaxy 15 Mpc away. 



Kervellaet al. [2005 1. 
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Younes, Porquet, Sabra, fc Reeves ||2011 ). However, the broad Ha line emission observed in 



2010 is dominated by the double-peaked broad line component that is attributed to the tidal 
disruption event. Thus, it is quite conceivable that ram pressure shock ionization produced 
by the interaction of the in-falling supergiant star with the ambient interstellar medium 
could provide an additional mechanical source of ionization that would help alleviate the 
ionizing deficit noted for NGC 4203. Note that ram pressure ionization is a very different 
mechanism for exciting the contrail than proposed by Rees [ (1988) which appeals to external 
flares generated by material as it is accreted onto the BH (e.g. Gezari et al. 2003). Ram 
pressure ionization is produced in-situ. Although ram pressure ionization was discussed 
previously in the context of the tail of ionized gas leaving the Galactic Center M supergiant 
IRS 7 (Serabyn, Lacy, & Achtermann 1991 Yusef-Zadeh & Melia 1992) there is a need to 
explicitly model the tidal disruption of a supergiant star as existing studies have tended to 
focus on solar type stars (e.g. Bogdanovic et al. 2004, and references therein). 



5. Conclusions 

Spectroscopic observations with the Hubble Space Telescope (HST) have revealed a time 
variable double-peaked broad Ha emission line profile in NGC 4203 which has been success- 
fully modeled as the combination of a time variable disk plus a steady state inflow. In this 
model the broad Ha line emission observed in 1999 is dominated by the inflow which is large, 
~ 2 pc in diameter, and likely sustained by stellar mass loss. If the gas density is high, > 
10^ cm~^, as suggested by the absence of similarly broad [O I] and [O HI] emission lines, 
then interpreting the broad Ha emission line observed in 1999 leads to a steady state inflow 
rate ~ 2 x 10~^ M0/yr which exceeds the requirement to explain the X-ray luminosity in 
terms of radiatively inefficient accretion by a factor of ~ 6. The time variable double-peaked 
component of the broad Ha emission line, which is particularly impressive in the 2010 spec- 
trum, has been modeled as an axis-symmetric disk, more appropriately described as a narrow 
ring, leading to the conclusion that it most likely represents the contrail of a supergiant star 
that is being tidally disrupted by the central BH. The AGN is apparently unable to sustain 
the ionization of the broad Ha emission line. The discrepancy is particularly acute in 2010 
unless the AGN had brightened by a factor of 4 since 2003 which seems unlikely. Ram 
pressure shock ionization produced by the interaction of the in-falling supergiant star with 
the ambient interstellar medium could provide an additional mechanical source of ionization 
that would help alleviate the ionizing deficit noted for NGC 4203. 

This research has made extensive use of the NASA Astrophysics Data System, the 
Atomic Line List, |http://www.pa.uky.edu/^peter/newpage/| and a variety STSDAS tasks. 
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Wavelength (A) 

Fig. 1. — Visual spectra of NGC 4203 as seen through the following gratings: Left panel: 
G430L. Right panel: G750M. Red line shows data obtained under PID 11571. Black lines 
for both panels show data obtained under PID 7361. 
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Fig. 2. — Broad Ha emission line in NGC 4203 observed in 1999 (Left panels) and 2010 
{Right panels). Top panels: The observed spectrum is shown in black and a model for the 
forbidden lines is shown in red (see also Table 2). Lower panels: The broad Ha emission 
line profile after the forbidden lines have been subtracted. The central wavelengths of the 
subtracted lines are indicated in red. The vertical black line corresponds to the observed 
(redshifted) central wavelength of the Ha line 
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Fig. 3. — Broad H/3 emission line in NGC 4203 observed in 1999. Top panel: The observed 
spectrum is shown in black and a model for the forbidden lines is shown in red (see also 
Table 2). Lower panel: The broad H/3 emission line profile after the forbidden lines have 
been subtracted. The central wavelengths of the subtracted lines are indicated in red. The 
vertical black line corresponds to the observed (redshifted) central wavelength of the H/3 line 
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Fig. 4. — Broad H7 emission line in NGC 4203 observed in 1999. Top panel: The observed 
spectrum is shown in black and a model for the forbidden line is shown in red (see also 
Table 2). Lower panel: The broad H7 emission line profile after the forbidden line has been 
subtracted. The central wavelength of the subtracted line is indicated in red. The vertical 
black line corresponds to the observed (redshifted) central wavelength of the H7 line 
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Fig. 5. — (Left Panel). Comparison of the normalized broad Ha, H/3 and H7 emission lines 
in NGC 4203 observed in 1999. 
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Fig. 6. — The broad Ha line observed in 1999 (Left panel) and 2010 {Right panel) modeled as 
a combination of a disk (blue line) and an inflow (red line). The thick black line represents the 
data. The sum of the model components and the residuals are plotted as thinner gray lines. 
These are just two of the many possible solutions permitted by the surface illustrated in 
Fig. 7. The model disks plotted here both have an inclination of 82° but different sizes and 
brightnesses (see Table 6) and represent the solutions at the tip and the tail of the arrow 
plotted in Fig. 7. 
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Fig. 7. — Contour plot of the reduced surface constraining the model disk parameters 
that best represent the NGC 4203 spectra observed in 2010 (red contours) and 1999 (blue 
contours). The thicker lines define the xLin contour. The thinner lines define the xLin + 2c" 
contour. There are a number of degenerate solutions for each set of observations. However, 
the red and blue Xmm contours do not overlap suggesting that there is no single disk model 
that can simultaneously explain both the 1999 and 2010 observations. The arrow represents 
the least size evolution solution (see Section 4.7 for details). 
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Table 1. NGC 4203 Spectral Datasets 



PID 


Observation Date 


Grating 


Spectral Range 


Slit 


Dispersion 


Plate Scale 


Integration Time 


Datasets 








A 


arc sec 


A/pixcl 


arc sec/pixel 


s 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


7361 


4-18-1999 


G750M 


6295 - 6867 


52 X 0.2 


0.56 


0.05 


900 


o4e010010 


7361 


4-18-1999 


G750M 


6295 - 6867 


52 X 0.2 


0.56 


0.05 


979 


o4e010020 


7361 


4-18-1999 


G750M 


6295 - 6867 


52 X 0.2 


0.56 


0.05 


900 


o4e010030 


7361 


4-18-1999 


G430L 


2900 - 5700 


52 X 0.2 


2.73 


0.05 


730 


o4e010040 


7361 


4-18-1999 


G430L 


2900 - 5700 


52 X 0.2 


2.73 


0.05 


900 


o4e010050 


11571 


3-2-2010 


G750M 


6295 - 6867 


52 X 0.1 


0.56 


0.05 


860 


obSiOlOlO 


11571 


3-2-2010 


G750M 


6295 - 6867 


52 X 0.1 


0.56 


0.05 


860 


ob3i01020 


11571 


3-2-2010 


G750M 


6295 - 6867 


52 X 0.1 


0.56 


0.05 


728 


ob3i01030 
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Table 2. Emission Line Parameters for the G750M Nuclear Spectrum Obtained 4-18-1999'' 



Line 


Central Wavelength'^ 


Flux 


FWHM 




A 


10-15 erg cm-2 3-1 


kms— 1 


(1) 


(2) 


(3) 


(4) 


[O I] 


6325 ± 1 


12.3 ± 2.0 


630 ± 131 


[O I] 


6390 ± 3 


3.9 ± 2.1 


649 ± 440 


[N II] 


6574 ± 3 


6.0 ± 3.2 


400 ± 147 


Ha (broad) 


6588 


168 ± 2 


1300 ± 100 


Ha (broad)<^ 


6588 


372 ± 4 


2667 ± 100 


Ha (narrow) 


6588 ± 3 


5.0 ± 2.0 


300 ± 85 


[Nil] 


6609 ± 1 


21.2 ± 2.5'= 


400 ± 28 


[S II] 


6742 ± 1 


3.5 ± 0.7 


283 ± 75 


[S II] 


6757 ± 1 


5.1 ± 0.8 


371 ± 75 



^Measured within a 0.2" x 0.35" aperture. Continuum subtracted but not 
corrected for dust extinction. 

''Observed wavelength 

^=1999 

'^2010, measured in a 0.1" x 0.35" aperture 

"The [N II] emission line flux is chosen so as to not over-subtract the broad 
Ho emission line profile. The [N II] line shape is modeled using two components; 
a logarithmic profile of width 400 km/s and a gaussian of width 130 km/s each 
contributing 90% and 10% of the emission line flux, respectively. 
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Table 3. Emission Line Parameters for the G430L Nuclear Spectrum Obtained 4-18-1999^ 



Line 


Central Wavelength'' 


Flux<= 


FWHM 




A 


10"^'* erg cm~^ 




(1) 


(2) 


(3) 


(4) 


[O II] 


3742 ± 1 


0.96 ± 0.02 


764 ± 24 


H7 (broad) 


4360 


0.8 ± 0.1 


1770 ± 200 


[O III]d 


4384 ± 1 


< 0.1 


570 


H/3 (broad) 


4882 


3.33 ± 0.04 


1263 ± 100 


[O III]<= 


4981 ± 1 


0.9 ± 0.1 


570 


[O III]" 


5028 ± 1 


2.2 ± 0.2 


570 ± 64 



"Table entries that do not include uncertainties are fixed parameters. 
''Observed wavelength 

■^Measured within a 0.2" x 0.35" aperture. Continuum subtracted but not 
corrected for dust extinction. 

''The [O III] emission line flux is chosen so as to not over-subtract the broad 
H7 emission line profile 

"Adopting a logarithmic profile shape. 
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Table 4. Inflow Plus Disk Model Fixed Parameters 



Parameter Value 



Disk outer to inner radius ratio, (,o/£.i 3.75 

Disk Broadening parameter, a 350 km s~^ 

Disk Emissivity index, q 3.0 

Inflow inner radius, rin„er 0.01 pc 

Inflow outer radius, router 1 pc 

Black hole mass, M, 6 X 10'' M0 
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Table 5. Parameters for Computing Xred 





Observation 


Date 


Parameter 


1999 


2010 


Summation performed over j points 


461 


552 


Velocity range of summation, km/s 


-5860 to 5860 


-7010 to 7010 


Degrees of freedom, v 


452 


543 


Minimum , 

^red 


6 


1 


Inflow/Disk Intensity ratio normalized at v=0 


10 


2.75 



= j data points - 8 parameters - 1 



